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ABSTRACT 

Recently Gatewood, Han & Black ( 2001| ) presented an analysis of the Hipparcos stellar 
positions and their own ground-based measurements of p CrB, suggesting an astromet- 
ric orbit of 1.5 milli-arc-sec with an extremely small orbital inclination of 0.5°. This 
indicates that the planet-candidate secondary might be a late M star. We used the 
Hippa rcos data of p CrB together with the individual radial velocities of Noyes et al. 
( 1997 ) to independently study the stellar orbit and to assess its statistical significance. 
Our analysis yielded the same astrometric orbit. However, a permutation test we per- 
formed on the Hipparcos measurements indicated that the statistical significance of 
the astrometric orbit is only 2cr. Therefore, we can expect about one out of 40 systems 
with a similar data set to show a false astrometric orbit. 
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1 INTRODUCTION 



One of the first stars around which a planet candidate was found (Noyes et al. 1997) was p Coronae Borealis (HR 5968 = 
HD 143761 = HIP 78459). The discovery was based on the detection of a small periodic radial- velocity modulation, with an 
amplitude of 67 m s~^. The corresponding minimum mass of the unseen companion was found to be 1 Jupiter mass (=A/j), 
its actual mass depending on the orbital inclination. 

Recently Gatewood, Han & Black (2001) presented an analysis of the Hipparcos intermediate data, together with new 
ground-based astrometry they obtained with the Multichannel Astrometric Photometer (Glatewood 1987). Their analysis 
suggests an orbital semi-major axis of 1.5 milli-arc-sec(=mas) with an extremely small orbital inclination - 0.5°, indicating 
that the secondary might be a late M star with a mass of 0.14 ± 0.05 M©. The extremely small orbital inclination makes 
this discovery intriguing, because the probability of a binary to have such a small angle, assuming an isotropic orientation in 



space, is extremely small - 4 x 10 . Even among 50 systems, the current number of planet candidates (Schneider 2000), the 



probability to find one system with such an extremely small inclination is only 0.002. Therefore, the statistical significance of 
the Gatewood, Han & Black (2001) analysis, which might have implications beyond the study of p CrB, should be estimated 
carefully. 



Recently, Pourbaix (MOl) studied the significance of the derived Hipparcos orbits of many planet candidates, including 
p CrB. Pourbaix fitted an orbit to the Hipparcos data and checked the improvement in the fit resulting from the additional 
parameters, using an F-test. He used the F-distribution to assess the final significance level of the derived orbit, concluding 
that the significance is somewhat high - 99 per cent. This use of the F-distribution assumes Gaussianity of the individual 
measurements. We avoid this assumption by using the permutation test, which belongs to the class of distribution-free tests 
(e.g. Good 1994) and thus is more robust against modeling problems of the measurement process. 

We present in this paper an independent analysis of only the Hipparcos astrometric data (i.e. without the MAP data), 
with an emphasis on the assessment of the significance of the detection. On one hand, we show that the Hipparcos data alone. 
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Table 1. Best-Fitting Orbital Solution for p CrB. 



p 


39.73 ± 0.11 


days 


K 


67.8 ± 2.8 


m s ^ 


e 


0.023 ± 0.047 




To 


2, 449, 060 ± 16 


HJD 




7° ± 129° 




7 


47.8 ±2.2 


m s"-*- 


ai sin i 


2.48 ±0.10 


10-*AU 


i 


0? 54+0- 24 






-0?13 




n 


118° ±52° 




ao 


1.49 ±0.46 


milli-arc-sec 


parallax 


56.70 ±0.63 


milli-arc-sec 



together with the radial-velocity measurements of Noyes et al. (1997), yield, indeed, a small astrometric orbit of 1.5 mas, 
implying a small inclination for p CrB. On the other hand, we show that the statistical significance of this finding is only 
about 2a. Therefore, we can expect about one system out of 40 to have shown a false astrometric orbit. Section 2 presents 
our analysis, while Section 3 discusses our finding. 



2 ANALYSIS 

2.1 Orbital Solution 

The present analysis used the 38 AFOE radial velocities of p CrB (Noyes et al. 1997; 1999) , together with the 42 astrometric 



measurements of Hipparcos (ESA 1997), which were analysed by the FAST and the NDAC consortia (van Leeuwen & Evans 
199^). The spectroscopic and astrometric solutions have in common the following elements: the period, P; the time of periastron 
passage, To; the eccentricity, e; the longitude of the periastron, ui. In addition, the spectroscopic elements include the radial- 
velocity amplitude, K, and the centre-of-mass radial velocity 7. The astrometric orbital elements include three additional 
elements - the angular semi-major axis of the photocentre, ao; the inclination, i; the longitude of the nodes, fi. In addition, 
the astrometric solution includes the five regular astrometric parameters - the parallax, the position (in right ascension and 
declination) and the proper motion (in right ascension and declination). All together we had a 14-parameter model to fit to 
the spectroscopic and astrometric data. Note that we have allowed small corrections to the values of the Hipparcos reference 
solution of the five regular astrometric parameters. 

These 14 elements are not all independent. From K, P and e we can derive the projected semi-major axis of the primary 
orbit - a\^phys x sini, in physical units. This element, together with the inclination i and the parallax, yields the angular 
semi-major axis of the primary, ai. Assuming the secondary contribution to the total light of the system is negligible, this is 
equal to the observed ao- 

The results of our fit, which are consistent with the previous orbital solutions (Noyes et al. 1997, 1999; Gatewood et 
al. 2001), are given in Table |l| The small semi-major axis implies a secondary mass of M2 = 0.125 ± 0.042 Mq, where we 
assumed 1 A/© for the primary ( [Noyes et al. 1997| ). 



2.2 Significance 



The semi-major axis of the derived astrometric orbit is 1.49 ± 0.46 mas. A very small semi-major axis could have been falsely 



"detected" even without any real astrometric motion, due to the scatter of the actual measurements (Halbwachs et al. 2000) 



To find out if this is the case here we performed a "permutation" test (e.g., Good 1994) by generating simulated data from 
the very same astrometric measurements of p CrB. If there is some evidence of an orbit in the measurements, it should be 
ruined by the permutation. However, if the derived orbit is spurious, some random permutations should be able to reproduce 
a similar effect. In a sense, we let the data "speak for themselves" and do not have to assume any specific distribution for the 
measurements. 



We used the IAD Hipparcos measurements ( ESA 1997 ) and per muted the a ctual timing of the observations, modifying the 
partial derivatives with respect to the five astrometric parameters ( ESA 1997 ) accordingly. We then analysed the permuted 
astrometric data together with the actual radial velocities, deriving a new false astrometric orbit. 

For most of the original measurements there are two stellar positions, one derived from the NDAC and the other from 
the FAST consortia. The two positions are, obviously, not independent, but have an assigned non- vanishing correlation. In 
our permutation we kept the pairing of the corresponding NDAC and FAST positions, while permuting the timings among 
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Figure 1. Histogram of the size of the falsely "detected" semi-major axes in the simulated permuted data. The size of the actually 
detected axis is marked by an arrow. 



the pairs. This turned out to be a nontrivial exercise, and caused us to falsely assign a higher significance to the astrometric 
detection in an early phase of this study. 

A histogram of the falsely detected semi-major axes, derived from the "simulated" permuted data, is presented in Figure 
^ We see that most of the artificially detected semi-major axes are at the range of 0.5-1.0 mas, with a small distribution tail 
beyond 1.5 mas. In 2,000 random permutations we got 46 astrometric solutions with semi-major axis larger than 1.5. This 
indicates that the detection significance is at the level of 0.977, which is about 2cr. 



3 DISCUSSION 

The analysis presented here shows that on average we expect about one out of 40 systems with a similar dataset to p CrB to 
show a false astrometric orbit with a = 1.5 mas. This implies that the statistical significance of the claimed astrometric orbit 
is not high enough at this point of the research. 



As explained in the introduction, Pourbaix (2001) derived that the significance of the astrometric solution is at the 99 
per cent level. Compared to his technique we use a distribution- free analysis and thus we feel that our assertion, of a 97.7 per 
cent level, is more justified. 

One way to find out the true nature of the secondary of p CrB is to try and observe the secondary directly. If the secondary 
is indeed a stellar obje ct, infrared spectroscop ic observations, especially when analysed with TODCOR - a two-dimensional 
correlation technique ( IZucker fc Mazeh 1994| ), should detect some trace of the faint secondary (e.g., Mazeh et al. 2000a). 
Together with M. Simon and L. Prato, we plan such observations in the near future. Obviously, more precise astrometric 
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The Planet Candidates 




Period [days] 

Figure 2. The stellar metallicity as a function of tlie planet-candidate orbital period for planets with periods shorter than 100 days. 
The point representing p CrB is circled. 



measurements, expected to come in the future from the FAME (Horner et al. 1999), SIM (NASA 1999), DIVA (Roser 199E 



and GAIA (Gilmore et al. 1996) missions, will enable us to derive the secondary mass much more accurately. 



We turn now to discuss the relatively low metallicity of p CrB. It was noted by many workers (Gonzalez 1997; Marcy 
& Butler 1998; Queloz et al. 2000; Gonzalez 2001; Butler et al. 2000) that most stars that were found to harbor planet 



candidates exhibit metallicities higher than the typical metallicity of the solar neighborhood. Queloz et al. (200C) and Butler 
et al. ( 200C| ) further pointed out that the host stars to the "51 peg like" planets are particularly metal rich. Mazeh & Zucker 
(2001) suggested a probable decrease of the stellar metallicity as a function of the orbital period of the planet candidates, a 
dependence which holds, according to their suggestion, up to about 100 days. To study their suggestion we plot in Figure ^ 
the metallicity as a function of the planet orbital period, for the stars that appear in The Extra Solar Planets Encyclopedia 



web site as of September 2000 (Schneider 2000). We follow the Mazeh & Zucker (2001) suggestion and therefore include in 



the figure only planet candidates with periods shorter than 100 days. 

Whenever available we have used metallicities derived from spectral analysis, mostly from the seminal work of Gonzalez 



( 2001 ). Mazeh et al. (20001) derived the metallici ty of HP 209458. Whenever such an analysis was not available, we have used 
the photometric metallicity derived by Gimenez (2000). The metallicities of the stars not considered by Gimenez were derived 
by us follo wing his prescript ion, b ased on the photometry of Hauck & Mermilliod (|l99^ ) and the calibrations calculated by 



Crawford (|l97q ) and Olsen (|198J. 

In this figure p CrB has one of the lowest metallicities. The analysis presented here, which does not support the conjecture 
about the stellar nature of the companion of p CrB, strengthens the suggestion that the metallicity might depend on the 
period. The figure draws attention to the other two low-metallicity planet candidates. The extreme one is HD 114762 (Latham 
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et al. 1989; Mazeh, Latham & Stefanik 1996) with a metallicity of about —0.60. The nature of this secondary is not clear 
because it has a minimum mass of about 10 M,j, and even a moderate inclination can turn it into a brown-dwarf or stellar 
secondary. However, The other one is HD 6434 ( [Udry et al. 2000| ), with a minimum mass of 0.5 Mj. This small minimum 
mass renders the conjecture that this is a real planet safe. Therefore, even if the secondary of HD 114762 turns out to be 
stellar, the dependence of the metallicity on the period might still be real. 
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